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ABSTRACT 

We present observations and magnetic field modeling of the large polar crown 
prominence that erupted on 2010 December 6. Combination of SDO/AIA and 
STEREO_Behind/EUVI allows us to see the fine structures of this prominence 
both at the limb and on the disk. We focus on the structures and dynamics of this 
prominence before the eruption. This prominence contains two parts: an active 
region part containing mainly horizontal threads, and a quiet Sun part containing 
mainly vertical threads. On the northern side of the prominence channel, both 
AIA and EUVI observe bright features which appear to be the lower legs of 
loops that go above then join in the filament. Filament materials are observed 
to frequently eject horizontally from the active region part to the quiet Sun part. 
This ejection results in the formation of a dense-column structure (concentration 
of dark vertical threads) near the border between the active region and the quiet 
Sun. Using the flux rope insertion method, we create non-linear force-free field 
models based on SDO/HMI line-of-sight magnetograms. A key feature of these 
models is that the flux rope has connections with the surroundings photosphere, 
so its axial flux varies along the filament path. The height and location of the dips 
of field lines in our models roughly replicate those of the observed prominence. 
Comparison between model and observations suggests that the bright features on 
the northern side of the channel are the lower legs of the field lines that turn into 
the flux rope. We suggest that plasma may be injected into the prominence along 
these field lines. Although the models fit the observations quiet well, there are 
also some interesting differences. For example, the models do not reproduce the 
observed vertical threads and cannot explain the formation of the dense-column 
structure. 
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1. INTRODUCTION 

Solar prominences are relatively cool structures embed ded in the million- degree corona 
( IHirayamal Il985l ; lLabrosse et al.l |2010| ; iMackay et al.l 120101 ) . In Ha when viewed above the 
solar limb, prominences appear as bright structures against the dark background, but when 
viewed as "filaments" on the solar disk they are darker than their surroundings. We will use 
the terms "filament" and "prominence" interchangeably in general. A filament is formed in a 
filament channel, which is defined as a region in the chromosphere surrounding a polarity in 



version line (PIL ) where the chromospheric Ha fibrils are aligned with the PIL ( ]Foukallll971 



Gaizauskasl Il998l ) . Filaments can be found inside activity nest consisting of multiple bipo- 
lar pairs of spots ("active region filaments"), at the border of active regions ("intermediate 
filaments"), and on the quiet Sun ("quiescent filaments"), including the polar crown. 

Filaments typically consist of three structural co mponents: a spine, barbs, and two 



extreme ends (jMartinlll998l ; iLin et al.ll2008l ; iLinl 1201 ll ). The spine runs horizontally along 



the top of the filament, although there may be sections along the filament where the spine 
is nearly invisible. The barbs protrude from the side of the filament and when observed 
closer to the limb, the barbs are seen to extend down from the spine to the chromosphere 
below. The ends, also called "legs", may be a collection of threads that appear to termi- 
nate at a single point or at multiple points. When a quiescent filament is viewed with high 



resolution on the solar disk, Ha observations indicate that each o: 



: these three structural com- 



ponents consist of thin thread-like structures ( ILin et al.ll2008 
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The magnetic field plays a primary role in filament fo rmation, stability, and eruption 



( jPriest Ill989l ; iTandberg-Hanssen Ill995l ; IMackay et al.ll2010l ). However, the magnetic struc 



ture of prominences is still not fully understood, with many observations and theoretical 
models differing on the exact nature of the magnetic field. Various models for prominence 
magnetic structure can be summarized as follows. In most of the models, filament plasmas 
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these models, the prominence plasmas are embedded in a large sheared arcade or a flux 
rope that lie horizontally above the PIL. Flux rope models may be split into two categories: 
weakly twisted flux rope and highly twisted flux rope. Weakly t wisted flux rop e mode ls are 
similar to sheared arcades with one key difference as argued by iMackay et al.l ( 120101 ) . For 
flux rope models the flux rope and overlying arcade are independent flux systems with a 
separatrix surface between them. In contrast, for a sheared arcade model only a single flux 
system exists. In addition to the above mentioned 3D global models, there are also two other 
dip models which focus on the loc al support of vert i cal threads in hedgerow prominence, i.e., 
"sagged horizontal fields model" (|Chae et al.ll2008l ; [Chad 12010) and "tangled magnetic field 



model" (Ivan Ballegooijen fc Cranmerl 120101 ) . In the earlier model, the vertical threads are 
stacks of plasma supported against gravity by the sagging of initially horizontal magnetic 
field lines. In the latter model, tangled fields exist in a vertical current sheet of small width 
that is confined by the vertical fields on either side of the sheet. Neither of these models 
describes the global 3D topology of the magnetic field supporting the fila ment. In contras t 
to models mentioned above, the "field aligned thread model" developed by lLin et al.l ( 120081 ); 
Linl ( 120111 ) does not contain dips. This model is a 3D empirical magnetic model based on 
high-resolution Ha observations and is based on the assumption that the observed fine scale 
structures are parallel to the magnetic field. In this model, the filament plasma is located on 
magnetic arches that are highly sheared in the direction along the PIL. Some field lines run 
along the entire length of the filament and outline the filament "spine" . Other shorter ones 
run partially along the spine, but spread out from it and connect down to minority polarity 
elements on either side of the PIL. These shorter structures represent the filament barbs. In 
this case, the plasma is supported by MHD waves. However, relatively high frequencies and 
wave amplitudes are required, and it is unclear why such waves would not lead to strong 
heating of the prominence plasma. 

The purpose of the present paper is to analyze observations of a quiescent prominence, 
and to develop 3D models of its magnetic structure. Filaments in active regions have been 
modeled extensively in the past. Several authors have developed non-line ar force-free field 
(NLF FF) models which are based on magnetic observations. For example, iRegnier &: Priest 
(120071 ) construct N LFFF models by "ex t rapolating" observed p hotospheric vector fields into 
the corona (also see lRegnier et al.ll2002l ; ICanou fc Amarill2010l ). van Ballegooijen (2004) de- 
veloped an alternative method for constructing NLFFF models of filament flux ropes. The 
method involves inserting a flux rope into a potential field based on an observed photospheric 
magnetogram, and then evolving the field in time to an equilibrium s tate u sing magneto 
frictional relaxation. This method has been used by Bobra et al.l (120081 ) and ISu et al. 



J2009al lbl. | 20111 ) ; lAsgari-Targhi fc van Balleeooiienl (120121) to study active reg ions with fil- 
aments, and by iSavcheva fc van Ballegooijen! (120091 ); ISavcheva et al.l (j2012al Jbl) to study the 
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evolution of soft X-ray sigmoid. The papers mentioned above suggest that the observed 
active region filaments can be well explained by a flux rope model. In this paper we apply 
this flux rope insertion method to a quiescent prominence for the first time. 



Observations 



2.1. Data Sets and Instruments 



A large polar crown prominence w as observed in the 



Atmospheric Imaging Assembly (AIA, iLemen et al. 



(SDO), as well as the STERECLB (Behind) /EUVI flWuelser et al 



period December 1-6 in 2010, by the 
2012|) aboard the Solar Dynamics Obse rvatory 
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Golub et al 
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Howard et al 
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Kano et al. 



2008 



2008) 



Synoptic observ ations by the X- Ra y Telescope (XRT, 
aboard Hinode ( iKosugi et al.l 120071 ) and Ha observations by the Kanzelhohe Solar Observa- 
tory (KSO) are also included in the study. The photosphe ric magnetic field information is 
provided by the Helioseismic and Magnetic Imager (HMI, ISchou et al.ll2012l ) aboard SDO 
and SOLIS (Synoptic Optical Long-term Investigations of the Sun). 

The large prominence erupted around 14:18 UT on 2010 December 6. This prominence 
eruption is associated with a CME, which has a linear speed of 538 km s _1 , according to the 
SOHO LASCO CME catalo g. D e tailed studi e s on the prominence eruption are presented 



by ISu fc van Ballegooijenl (j2012bl ): iThompsonl (120121 ) . In the current paper, we present the 
structure and dynamics of the prominence before and after the eruption. 



2.2. Structure and Dynamics of the Prominence and Surroundings 

2.2.1. Prominence Structure and Dynamics before the eruption 

Figures [TH2] shows SDO/AIA observations of the evolution of the prominence within 
six days (December 1-6) prior to the eruption. The images on the left, middle, and right 
columns are taken at 304 A, 171 A, and 193 A, respectively. Note that artificial color is used 
for the on-line only figures. Red color is used for the images taken at 304 A by both AIA 
and EUVI, while 171 A images by AIA are yellow. Images at 193 A by AIA are dark yellow, 
while EUVI images at 195 A are green. Moreover, both 193 A and 195 A provide similar 
spectral information. On December 1, the prominence is observed at the east limb and 
contains mainly horizontal threads shown as emission at 304 A (white arrows in Figure [JjV). 



1 http://cdaw. gsfc.nasa.gov/CMEJist/ 
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A dark feature is also seen in absorption at all three wavelengths (black arrows in Figures 
[]Ji4Tb)- On December 2 and 3, the first part of the prominence rotated onto the disk, and the 
prominence at the limb becomes more complex. One can see some dark vertical threads with 
overlying bright horizontal threads as shown in Figures Id and lg. Images at 304 A shows 
the vertical threads either in emission or absorption. The dark vertical threads are better 
seen at 171 A and 193 A. Sometimes, bright emission in the promin ence is also vi s ible a t 



171 A (e.g., Fi gure[Th). This bright emission in 171 A is analyzed by IParenti et al.l ( 120121 ). 
Figure [2] shows the prominence from December 4 to December 6. The filament on the disk 
shown at 171 A and 193 A appears to be much narrower than the one shown at 304 A. 
SDO/HMI observations (black and white contours in Figure [2ji) show that the underlying 
magnetic fields are positive on the southern side and negative on the northern side of the 
filament. 

The AIA observations suggest that this prominence can be divided into three parts: a 
region with mostly horizontal threads (Part I), a region with horizontal threads overlying 
the vertical threads (Part II), and a region with mostly vertical threads (Part III). This 
prominence contains several different features which are marked as numbers in Figure [2g. 
Part I is located in the active region, which is on the northwestern side of the black dashed 
line in Figure |2g. The horizontal threads (marked as "1" in Figure |2g) in Part I are best 
viewed at the limb in Figured^. Part II is located between the black and white dashed lines 
as shown in Figure [2g. The limb view of this region is shown in the images from December 
2 to December 6 in Figures dH2J Feature "2" is a dense column composed of dark vertical 
threads which appear to be perpendicular to the overlying horizontal threads. The limb 
view of Part II displays an arch-like structure. The two feet of this arch (marked as "4" and 
"6") contain mostly quasi vertical threads, while the top of the arch is composed of mostly 
horizontal threads (feature "3"). Note that the threads in feature 4 are not vertical, and 
the prominence spine appears to tilt towards the South. Feature 5 refers to the hole located 
below the top of the arch. Part III is located on the northern side of the white dashed 
line. This part is composed of mostly vertical threads (feature "8") overlying a bubble-like 
structure (feature "7"). Feature "9" refers to a horn-like structure located above the vertical 
threads. All prominence structures are clearly visible at 304 A, while only parts of the 
prominence are seen as absorption in 171 A and 193 A, which best reveal the fine structure 
of the vertical threads. Note that some features are also visible as emission in 171 A. 

The aforementioned AIA observations reveal the filament while it rotating from the east 
limb to the disk. STEREO_B views this filament as it rotating from the disk to the west limb. 
Figures [SHU show the filament from December 1 to December 6 at 304 A and 195 A observed 
by EUVI aboard STEREO_B. On December 1, the filament displays a dark long and narrow 
structure. The filament on the east appears to grow wider and thicker with time. The black 
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dashed line in Figure H] divides the filament into two parts. The western part of the filament 
(active region part) contains mainly horizontal threads, while the eastern part (quiet Sun 
part) contains mainly quasi-vertical threads. Unlike AIA, EUVI observations do not reveal 
the horizontal threads located on top of the quasi- vertical threads. A comparison of Figures 
Hfc andHk shows that the active region part of the filament appears to rise up on December 6. 

Figure |5] shows the filament observed at Ha by KSO within four different days prior 
the eruption. Ha observations show similar filament structures as those found at the EUV 
channels by AIA. A comparison of Figures \5jp and EH shows that most part of the active 
region filament (marked as black arrows) disappeared before the eruption on December 6. 
This is consistent with the rise of the active region filament as observed by EUVI. 

The prominence before the eruption is very dynamic. Horizontal counter-streaming 
along the top of the arch as shown in Figure [2g is clearly seen in the AIA 304 A observations 
(see video 1). Downflows are also observed in part of the vertical threads. Part of the vertical 
threads are just waving around without clear flows as shown in the AIA 193 A observations. 
EUVI observations (video 2) show that the filament material is ejected from the active 
region part to the quiet Sun part. This ejection leads to the horizontal oscillations of the 
vertical threads and the formation of a dark tree-like structure (marked using white arrows 
in Figure H|) near the border between the active region and the quiet Sun. 

2.2.2. Dense Column Structure 

Figure E] shows comparisons of the AIA and EUVI observations of the filament on 
December 6. The field of view of the AIA and EUVI images are shown as black boxes 
in Figure [2} and Figure 0J, respectively. Figures [TH2] and Figures |3]|U suggest that the 
appearance of the filament is very different in AIA and EUVI observations. Therefore, to 
compare different filament features observed by the two instruments, we first identify a series 
of bright-point like features, which are marked as numbers in Figure El A careful comparison 
shows that the dense column observed by AIA corresponds to the tree-like structure observed 
by EUVI. This dense column is formed by the strong accumulation of quasi-vertical threads 
(marked as white arrows in all images). 

2.2.3. The Structure around the Prominence 

The prominence is located in a very narrow filament channel, which is best observed 
at 195 A by EUVI as shown in Figure |3j Straight faint features (traced using white lines 
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in Figure |3]F) are observed on the two sides of the channel. These features may refer to the 
lower legs of the overlying coronal arcades. Near the quiet Sun part, one can see some bright 
curved features (traced using black lines in Figure (3F), which are located on the northern 
side of the filament channel. No clear counterparts of these features can be identifi ed on the 



southe rn side of the channel. This emission asymmetry that was first identified by ISu et al. 



( 120101 ) is also clearly visible in AIA observations, see Figures |2} and|6j However, these bright 
curved features appear to be straight in the AIA view. To have a better understanding on 
this asymmetry, we trace these features back when they are observed at the east limb by 
AIA. Figure [7] shows AIA 171 A images on December 4-6. Each image is the average of 
30 images taken within 6 minutes. All images are contrast enhanced using a radial filter 
technique developed by S. Cranmer (2010, private communication). The top images show 
bright loops which start from the northern side of the filament channel (white arrows), arch 
over the filament, and finally disappear behind the filament or join in the filament. Evolution 
of this type of loops from December 3 to December 6 are shown in the online video (video 
3). The existence of these loops is consistent with the emission asymmetry in the s ense that 



they d o not have clear footpoints on the southern side of the filament channel. ISu et al. 



( 120101 ) interpreted these bright curved features as magnetic field lines that turn into the flux 
rope. 

The bottom images in Figure [7] show highly sheared bright loops (white arrows) that go 
above the prominence. However, it is not clear where these loops end. Figure [7H also shows 
several bright horn-like features (black arrows) located right above the vertical prominence 
threads. This observation appears to suggest that some part of the flux rope lies above the 
prominence vertical threads. 

Figure [8] shows images of the prominence at different channels observed by SDO/AIA 
and Hinode/XRT on 2010 December 4. A clear cavity surrounding the prominence can be 
seen in the two XRT images. The AIA image at 335 A which is an average over 4 minutes 
shows similar structure as observed by XRT. No clear cavity can be identified in the other 
three EUV channels, i.e., 304 A, 171 A, and 193 A. The AIA image at 335 A (Figure EH) 
shows that the dark filament lies on the southern edge of the cavity rather than at the 
center. We also can see that the bright coronal arcades surrounding the cavity appear to 
be asymmetric, i.e., the bright structure on the northern side is brighter than that on the 
southern side. This asymmetry is also clearly visible in the two XRT images (Figures [8k- 

EJ). 



Martin fc McAllisterl (119961 ) found that there is a one-to-one correspondence between 



the dextral/sinistral orientation of the filament channel and the skew of the overlying coro- 
nal arcade. This implies that the axial magnetic field within filaments points in the same 
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direction as that of the surrounding corona. As viewed from the positive polarity side, for a 
sinistral filament the axial component of th e field in the arcade points to the left, and vice 
versa (i.e., Figure 1 in McAllister et al.ll2002l ). Figure |H] shows the post-eruption arcade above 
the erupted filament observed at 335 A by AIA on 2010 December 7. As shown in Figure [2h, 
the positive polarity (white patches) is located on the southern side of the filament channel. 
This orientation of the arcade suggests that this filament is a sinistral filament, according to 
the aforementioned definition. 



2.2.4- Prominence Structure after the Eruption 

The prominence eruption on December 6 is a partial eruption, which means that part of 
the filament survived the eruption. Figure ITUl shows the east end of the quiescent prominence 
after the eruption. SDO observes this prominence at the east limb (left column), while 
STEREO_B views it on the disk (right column). The side view of this prominence by AIA 
suggests that it contains mainly vertical threads. Images at 304 A shows the vertical threads 
either in emission or absorption. Dark quasi-vertical threads are seen in absorption at both 
171 A and 193 A. In 171 A, one can also see a bright horizontal structure located on top of the 
dark vertical threads. The filament displays a thin dark structure, when viewing it from the 
top as shown in the STERECLB observations. The width of the filament at 304 A (Figure 
fTUb) is larger than that at the other two EUV channels (Figures fTOH^TOb) . Without the 
different point of view from SDO, it is very difficult for us to imagine that this thin dark 
filament contains vertical threads. It looks very similar to the active region filament which 
contains mainly horizontal threads. This thin dark filament structure may be the horizontal 
filament threads located on top of the vertical threads. Another possibility is that it is just 
the accumulation of the lined up dark vertical threads, and there are no horizontal threads. 
The spatial resolution of EUVI does not allow us to distinguish whether there are overlying 
horizontal threads or not. 



Models for the Observed Prominence 



It has been suggested that Quiescent Prominences (QPs) are locat ed in helical flux ropes 



that lie horizontally aboy e filament channels on the quiet Sun (e.g., iRust fc Kumar! 11994 



Low &: Hundhausenl 1995 ). The prominence plasma is assumed to be located at dips in the 
helical field lines. In this section we construct three-dimensional (3D) flux rope models for 
the observed prominence, then we compare the models with AIA and EUVI observations. 
The coronal arcade overlying a QP flux rope is thought to play an important role in its 
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equilibrium and stability. To have stable equilibrium, the magnetic tension of the overlying 
arcade must balance the magnetic pressure of the flux rope. If the overlying fields are too 
weak, the flux rope will rise and expand into the heliosphere. Therefore, an important open 
question is whether QP flux ropes can exist in stable equilibrium with their surroundings. 
The overlying fields are anchored in the quiet photosphere, which often has mixed polarity 
with only a small excess of one polarity over the other. It is unclear how effective such mixed- 
polarity fields can be in anchoring QP flux ropes. The purpose of the present modeling is to 
determine whether the observed QP can be explained in terms of a stable flux rope model. 



3.1. Methodology 

In this subsection we describe how the models are constructed. In the last few years, we 
have developed various tools for modeling non-potential magnetic fields in the solar corona, 
called the Coronal Modeling System (CMS). Its main purpose is to construct NLFFF models 
of the coronal magnetic field. The models are constructed by inserting a magnetic flux rope 
into a potential-field model of the region and then applying magneto-frictional relaxation 
(MFR). This "flux rope insertion" method i s quite flexible and provides information about 



the stability of the resulting magnetic fields (ISu et al.ll201l[ ). In the present study we apply 



this flux rope insertion method to a quiescent prominence for the first time. 

The observed QP is more than a solar radius in length, and the model needs to in- 
clude the overlying coronal arcade, so a large simulation domain with spherical geometry 
is required. However, we also need to resolve the mixed-polarity fields on the quiet Sun, 
which requires high spatial resolution (a few Mm or better). To achieve these goals, the 
CMS code uses variable grid spacing: at low heights the grid spacing is about 0.002 cos A 
R , where A is the latitude, but at larger heights the cell size is increased by powers of 2. 
The high-resolution (HIRES) computational domain extends about 117° in longitude, 48° 
in latitude, and up to 2.41 R from Sun center. The magnetic field B(r) in this domain is 
described in terms of vector potentials, B = V x A. 

The lower boundary condition for the HIRES region is derived from line-of-sight (LOS) 
photospheric magnetograms obtained with the SDO/HMI. Since the prominence is observed 
near the east limb, we have to use magnetograms that are taken several days after the promi- 
nence eruption on December 6. We combine four magnetograms taken on 2010 December 
8, 9, 10 and 11 (each at 14:00 UT) to construct a high-resolution map of the radial compo- 
nent B r of magnetic field as function of longitude and latitude at the lower boundary of the 
HIRES region. We also use a SOLIS synoptic map of B r to compute a low-resolution global 
potential field, which provides the side boundary conditions for the HIRES domain, and also 
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allows us to trace field lines that pass through the side boundaries of the HIRES region. 

The HIRES magnetic map is shown in Figure [TTa . Note that outside the active region 
the field has mixed polarity with dominantly positive polarity on the south side of the PIL and 
negative polarity on the north side. Based on this map alone it is difficult to recognize exactly 
where the PIL is located. Therefore, we used co-aligned AIA images and HMI magnetograms 
to locate the base of the prominence/filament on the magnetic map (blue curve in Figure 
[TTT) . This curve is the path along which the flux rope will be inserted into the model. At the 
two ends of the path (blue circles) the flux rope is anchored in the photosphere. Figure [TTb 
shows the zoomed-in view of magnetic map shown in Figure [TTh (red and green contours) 
overlaid on the HMI LOS magnetogram taken at 14:00 UT on December 6 (black and white 
image). This figure shows that there is no significant differences between the magnetic fields 
taken before the eruption (on December 6) and after the eruption (December 10 and 11). 
One may notice some shift in individual magnetic patches, but no significant flux emergence 
and cancellation can be identified. Moreover, the polarity inversion line appears to remain 
the same. It suggests that using HMI data several days after the eruption to model the 
prominence just before the eruption is reasonable. 



The methodology of flux rope insertion has been described elsewhere (e.g., iBobra et al. 



20081 ; ISu et al.ll2009bl ). so only a brief summary will be given. First, the potential field is 
computed from the HIRES and global magnetic maps. Then by appropriate modifications 
of the vector potentials a "cavity" is created above the selected path, and a thin flux bundle 
(representing the axial flux of the flux rope) is inserted into the cavity. Circular loops are 
added around the flux bundle to represent the poloidal flux of the flux rope. Then MFR is 
applied to drive the magnetic field toward a force-free state. MFR is basically an evolution 
of the magnetic field according to the magn etic induction equation with the velocity take n 



to be proportional to the Lorentz force (e.g., lYang et al.lll986t Ivan Ballegooijen et al.ll2000l ). 



The magnetic diffusion is kept as small as possible, so that the magnetic topology of the 
field is approximately conserved (some diffusion is needed to prevent numerical artifacts). 
If the axial and/or poloidal fluxes of the flux rope are too large, the model will not reach 
an equilibrium state but the field will keep expanding and moving to large heights, similar 
to what happens on the real Sun during an eruption. In this paper we are interested only 
in stable configurations, so we seek values of the model parameters that produce stable 
NLFFFs. In previous work we usually iterated the induction equation about 30,000 times, 
but we found that for QP flux ropes this is not enough to reach an equilibrium state. For 
the models presented here we use 90,000 iterations. 

For the models discussed in this paper, two magnetic sources with fluxes of ±2 x 10 20 
Mx were added to the 5 r -map at the two ends of the filament path in order to enhance 
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the stability of the flux rope. The axial flux near the ends of the (blue) path equals the 
flux of the sources ($ ax i = 2 x 10 20 Mx). The axial field points to the left (east) as seen 
by an observer looking at the QP from the south polar region, which has positive polarity. 
Therefore, the flux rope has sinistral orientation, consistent with the observed skew of the 
post-erup tion arcades (see Fi gure [9]) and with the hemispheric pattern of chirality of filament 
channels (IMartin et al.l 119941 ) . 

The observations clearly show a filament within the active region, and the STEREO_B 
observations provide constraints on its height. Preliminary modeling showed that in order 
to reproduce the observed active region filament the axial flux of the flux rope must exceed a 
certain minimum value: $ ax i > 6 x 10 20 Mx. However, such values were found to be too high 
for the eastern part of the filament overlying the quiet Sun (longitude < in Figure [TT]). 
Axial fluxes in excess of 4 x 10 20 Mx lead to the eruption of the eastern part as the surrounding 
quiet-Sun fields are too weak to anchor the flux rope to the photosphere. Therefore, we must 
consider models in which the axial flux $ ax ; varies with position along the flux rope. This is 
achieved by adding magnetic flux to (or removing flux from) the flux rope at certain places 
along the filament path. These places are indicated by the yellow "barbs" in Figure [TTJ In 
fact, real filaments hav e barbs th at may indeed be sites where magnetic flux enters or leaves 



the main filament (e.g.. lLinll201ll ). In Figure (TTJ the circle at the end of a yellow line segment 



indicates the source (or sink) of the added flux. As we follow the (blue) filament path from 
right to left in the figure, the two barbs near longitude +40° each add 3 x 10 20 Mx to the 
flux rope, and the three barbs near longitude +5° each subtract 2 x 10 20 Mx. Therefore, 
the axial flux in the active region (in between the two sets of barbs) is 8 x 10 20 Mx, and the 
axial flux on the quiet Sun (outside the two sets of barbs) is 2 x 10 20 Mx. The barbs near 
+5° are put on the south side of the PIL in order to force the turned-out field lines to arch 
over the flux rope on their way to the negative polarities on the north side. We find that 
this initial setup leads to stable models that are in rough agreement with the observations. 



3.2. Modeling Results and Comparison with Observations 

In this paper we focus on models representing possible magnetic configurations just prior 
to the eruption that started at 14:20 UT on 2010 December 6 (as discussed in section [3~T| the 
models are actually based on magnetograms taken several days later). Two values for the 
poloidal flux F po i of the flux rope are considered: 1 x 10 10 Mx cm -1 (Model 1), and 2 x 10 10 
Mx cm -1 (Model 2). Figure [121 shows results from the NLFFF models and a comparison 
with AIA observation. The top panels show the two models, projected onto the plane of the 
sky as seen from SDO. The bottom panels show co-aligned images of the Sun taken with 
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AIA in the 193 A and 304 A passbands. The colored curves in the top panels are magnetic 
field lines traced through the model, and the blue features indicate the locations of dips in 
the field lines, i.e., sites where the field lines are locally horizontal and curved upward (light 
blue for low- lying dips, darker blue at larger heights). The red and green contours show 
the magnetic flux distribution on the photosphere. Note that the field lines in Figure [i~2b 
are more twisted than those in Figure [T2k . which is due to the difference in the poloidal 
flux of the inserted flux rope. Both models reproduce more or less the observed prominence 
height (about 70 Mm) as seen in the AIA 304 A image (see Figure [T2"H). This means that 
the modeled flux rope is sufficiently thick, and the axis of the flux rope lies at sufficiently 
large height to provide a reasonable model for the observed prominence everywhere along 
its length. In the remainder of this paper we only consider the model with poloidal flux 
F pol = 10 10 Mx cm" 1 , i.e, Model 1. 

Figure [13] shows a zoomed-in view of the middle part of Figure [T2k . At this resolution 
the blue features resolve into small blue dots that indicate the locations of field-line dips. The 
dots lie on surfaces where the radial component of magnetic field vanishes, B r (r,6,<p) = 0, 
where r, 9 and are spherical coordinates. Note that the dots in Figure [T3] make patterns 
of horizontal and quasi-vertical lines. These patterns should not be taken too seriously 
because they are an artifact of the way the dips are chosen. Specifically, we locate the 
intersections of the surface B r = with the "edges" of the 3D grid on which the magnetic 
field B is defined (actually, the dots are short line segments indicating the direction of the 
horizontal field at the dips, but this is not visible in Figure [T3"j) . Therefore, the patterns 
of horizontal and vertical lines are created by interference with the 3D grid. However, we 
also see "ridges" in Figure 13 where the density of dots is higher than elsewhere. This 
indicates that the surface of dips is strongly warped, and in some places our line-of-sight is 
tangential to this warped surface. The warping of the B r = surface is due to the un even 



distribution of magnetic sources on the photosphere (e.g.. lAulanier fc Demoulinl 120031 ). It 



has been suggest ed that the quasi-vertical threads in quiesc ent prominences may be due to 



such ridges (e.g., iHeinzel &: Anzerl 120061 ; iDudfk et al.ll2008l ). However, we find that in our 



model the ridges are not always vertical, and the number of ridges per unit length along the 
prominence is far too small for the ridges to be identified with the observed vertical threads, 
which usually are densely packed. Therefore, the distribution of field-line dips in our NLFFF 
model does not reproduce the observed quasi-vertical threads. 

Figure [T41 shows a comparison of Model 1 with AIA and EUVI observations of the coronal 
arcade surrounding the prominence. The top panels show the same set of field lines seen from 
two different viewing angles: (a) from SDO, and (b) from STERECLB. The bottom panels 
show images in AIA 193 A and EUVI 195 A, which are dominated by emission in Fe XII 
195 A. These images show thin threads and fans that presumably are aligned with the local 
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magnetic field. On the quiet Sun the plasma density drops rapidly with height due to the 
gravitational stratification of plasma along the magnetic field lines. The Fe XII ion is formed 
at temperatures of about 1.5 MK, and the density scale height for such plasma in the low 
corona is about 75 Mm. The emissivity in Fe XII 195 A is proportional to density squared, 
which decreases even more rapidly with height. Therefore, we expect that in observations 
on the solar disk the Fe XII 195 A emission will be dominated by plasma in the lower parts 
of coronal loops. Accordingly, in Figures ir^b-TUiH we overplot the lower parts of the field 
lines as blue line segments shown in IT4"a-[T4"b; only heights less than 0.07 R Q (49 Mm) are 
shown. Close inspection of these images shows that the orientation of the line segments 
is roughly in agreement with the directions of the thread-like structures seen by AIA and 
EUVI. In particular, the model reproduces the bright features on the northern side of the 
PIL. Some of these threads are located on field lines that turn into the flux rope as we follow 
its path from east to west. This comparison between model and observations indicates that 
Model 1 reproduces the surrounding magnetic fields quite well. This model also confirms our 
previous interpretation of the emission asymmetry on the two sides of the quiescent filament 
channels. Those bright curved features on the north sid e of the PIL re present the lower part 
of the magnetic field lines that turn into the flux rope ( jSu et al.ll2010l ). 



Figure [TBI shows the magnetic field strengths at dips in the field lines, plotted as function 
of the height z of the dips above the photosphere. The prominence plasma is assumed to be 
located at such dips. All positions along the filament channel are included in this plot, except 
the two ends of the flux rope where the rope is anchored in the photosphere. The heights 
of the dips range up to 70 Mm, and the field strengths range from 40 G in the active region 
to 4.5 G in the quiescent part of the filament. The values of field strengt h are consistent 
with measurements of prominence magnetic fields using the Hanle effec t (e.g..lBommier et al. 
1994h . and with results from constant-a magnetohydrostatic modeling (lAulamer fe Demoulinl 
200 3|). The plasma pressure in quiescent pro minences lies in the range 0.1 - 2.0 dyne cm -2 
( llensen fe Wiiklll990t iHeinzel fe Anzerl 120061 ). so a magnetic field strength of 4.5 G corre- 
sponds to a plasma in the range 0.12 - 2.5, where (3 is the ratio of plasma- and magnetic 



pressures. IHeinzel fe Anzerl ( 120061 ) have shown that when /3 > 1 the weight of the promi- 
nence plasma will significantly distort the shape of the magnetic field lines; such effects are 
not taken into account in the present models. From the above estimates for j3 it is clear that 
gravity may significantl y deepen the fiel d -line dips in the quiescent prominence compared 
to our NLFFF model. IHeinzel fe Anzerl (120061 ) have developed 2D models of prominence 
threads supported by gravity-induced dips in the field lines. In the active region the effects 
of gravity are expected to be minor (/3 < 0.03). 

There are interesting differences between our models and observations. First, as men- 
tioned above the AIA observations show thin vertical threads in the prominence near the 
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limb, but our models do not have any such features in either the shapes of the magnetic 
field lines or the distribution of field-line dips. This is not surprising because our NLFFF 
models do not include the effects of p lasma pressure and gravity, which likely are responsi- 



ble for producing the vertical threads ( IKippenhahn fc Schluterlll957l ; iHeinzel &: Anzerl 12001 



3Sp 

m 



Low fc Petrid 120051 ; iHillier et al.ll2012l ). Second, the model cannot explain the formation of 



the dense-column structure. 

The PIL is shown as a blue curve in Figure [TTJ and as a red curve in Figure [T2"l i. Note 
that the filament spine is displaced to the South relative to the PIL in Figure IT2"l i. whereas 
the base of the filament appears to be located near the PIL. Therefore, the prominence does 
not lie vertically above the PIL, but it tilted to the South. This tilt is not reproduced in 
either Model 1 or Model 2. In the models the flux rope lies directly above the PIL, which is 
where the flux rope is inserted into the model. We create another model (Model 3) in order 
to see whether we can reproduce the tilt of the prominence by changing the path where the 
flux rope is inserted. Figure [16] shows top view of the field line dips (blue features) from 
Model 1 and Model 3. Model 3 has the same parameters as Model 1 except that the filament 
path (red curve) is displaced to the south in the longitude range between —35° and —5°. 
The path in Model 3 is along the filament spine rather than the filament base as used in the 
other two models (see Figure H2"H) ■ Figure [16b suggests that the field line dips are pushed 
back to the PIL (filament base) after 90000-iteration MFR, though the flux rope is originally 
inserted along the filament spine (red curve). Therefore, the tilt of the filament cannot be 
better reproduced by changing the path of the flux rope. This absence of tilt is also not due 
to any lack of convergence of the MFR; there simply is no force in the model to push the 
flux rope sideways. Therefore, the observed tilt of the QP is somewhat of a mystery. 

Figure [T7] shows comparison of the observed cavity and magnetic field lines from Model 
1. The AIA image in Figure [TTa shows that the cavity is tall and narrow rather than a half 
circle-like shape. The size and shape of the model flux rope is represented by the inner most 
dark blue line shown in Figure ITTb . A comparison of these two figures suggests that the 
observed cavity is taller than the modeled flux rope. The observations show that the bright 
features surrounding the cavity are nearly vertical with respect to the solar limb. However, 
the lower part of the modeled field lines on the south side of the cavity are much more 
inclined in comparison to the observations. 

The two models presented in Figure [12] lie at the edge of the stable regime. Indeed, an 
eruption does occur starting at 14:20 UT, shortly after the time for which the models are 
constructed. However, the eruption involves only the western part of the QP; the eastern 
part remains behind and is apparently quite stable (see Figure [10]). This partial eruption is 
difficult to understand on the basis of our models: the flux rope extends along the full length 
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of the QP, and therefore we expect that the eruption should remove the entire flux rope. 
This leads us to question the validity of the flux rope model for describing the magnetic 
structure of this QP. 



4. Discussion and Conclusions 

A large polar crown prominence partially erupted on 2010 December 6. To understand 
the magnetic support of the prominence, we focus on the structure and dynamics of the 
prominence before and after the eruption. The 6-day SDO/AIA observations of the filament 
near the east limb before the eruption suggest that this filament can be divided into 3 
parts. The active region part near the west end contains mainly horizontal threads. The 
middle part near the quiet Sun contains vertical threads with overlying horizontal threads. 
While the east end of the prominence that survived the eruption appears to be composed 
of mainly vertical threads. The post-eruption limb observations by AIA suggest that the 
prominence left behind contains mainly vertical threads, while a regular thin dark filament 
structure appears in the on-disk observations by EUVI. This thin dark filament structure 
may be the horizontal filament threads located on top of the vertical threads. Another 
possibility is that it is just the accumulation of the lined up dark vertical threads, and there 
are no hori zontal threads. Such st ructure, consisting of horizontally aligned threads, was 



such threads was discussed by iHeinzel fc Anzerl (120061 ) . 



reported by lSchmieder et al.l (120041) for another filame nt. Corresponding Ha appearance of 



STEREO/EUVI observes straight and faint features representing the lower legs of the 
overlying coronal arcades on the two sides of the filament channel. On the quiet Sun part, 
bright features with no clear counterparts are also identified on the northern side of the 
channel. This emission asymmetry is consistent with the existence of loops that originates 
from the northern side of the channel, then arch over the filament and disappear or join in 
the filament as shown in earlier limb observations at 171 A by AIA. A clear cavity corre- 
sponding to the filament channel is observed on December 4 at 335 A by AIA and XRT. The 
observations suggests that the filament is located on the southern edge of the cavity rather 
than at the center. The bright features (lower legs of the coronal arcades) on the north side 
of the cavity are brighter than those on the south side. The skew of the post-eruption arcade 
suggests that this filament is sinistral. 

The prominence before the eruption is very dynamic. AIA observes horizontal counter- 
streaming overlying the vertical threads. STEREO_B/EUVI observes filament material ejec- 
tion from the active region part to the quiet Sun part. This ejection leads to strong oscilla- 
tions of the vertical threads and the formation of a dense column near the border between 
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the active region and quiet Sun. The dense column refers to strong accumulation of dark 
vertical filament threads. 



Using the flux rope insertion method developed by Ivan Ballegooijenl (120041 ) we construct 
two non-linear force-free field models. The poloidal flux (1 x 10 20 Mx cm -1 ) in Model 1 is 
half of that in Model 2, so the flux rope in Model 1 is less twisted. In both models, the 
axial fluxes are 8 x 10 20 Mx and 2 x 10 20 Mx in the active region and quiet Sun, respectively. 
Since the prominence is observed at the east limb, the models are constructed based on 
the combined LOS photospheric magnetogram observed by SDO/HMI several day after the 
prominence eruption. Both models lie at the edge of the stable regime. 

Our models match some of the observed features quiet well. For example, the height 
and location of the field line dips in our models appear to match those of the prominence 
threads, though the narrow vertical thread-like structure is not reproduced in our model. 
The lower part of the magnetic field lines from our models can replicate the bright features 
on the two side of the filament channel quite well. Our model also confirms our previous 
i nterpretation of the emission asymmetry on the two sides of the channel on the quiet Sun 



(ISu et al.ll2010l ). This emission asymmetry is due to the fact that some of the bright features 
on the northern side (the brighter side) of the channel represent the field lines that turn into 
the flux rope. Therefore, no counterparts of these features can be found on the southern 
side of the channel. We suggest that plasma may be injected into the prominence along the 
aforementioned field lines that turn into the flux rope, since these field lines start from the 
photosphere. This is consistent with the thermal nonequilibrium me chanism for prominence 



formation as suggested by Antiochos, Karp en, and colleagues (e.g., lAntiochos fc Klimchuk 
1991c iKarpen et al.ll2005l ; iLuna et al.ll2012l ). In this model, localized heating above the flux 
tube footpoints produces evaporation of the chromospheric plasma, which condenses in the 
coronal part of the flux tube. 
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Fig. 1. — SDO/AIA observations of the quiescent prominence at 14:00 UT on 2010 December 
1 (top row), December 2 (middle row), and December 3 (bottom row). (A color version of 
this figure is available in the online journal.) 
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Fig. 2. — SDO/AIA observations of the quiescent prominence at 14:00 UT on 2010 December 
3 (top row), December 4 (middle row), and December 5 (bottom row). The white and black 
contours refer to the positive and negative polarities as observed by SDO/HMI at 14:00 UT 
on 2010 December 6. (A color version of this figure is available in the online journal.) AIA 
observations of the prominence at 193 A and 304 A within 24 hours before the eruption are 
also available (video 1) in the electronic edition of the Astrophysical Journal. The field of 
view of the online video contains only the prominence near the limb. The images are also 
rotated for 120 degree in the clockwise direction. 
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Fig. 3. — STEREO_BB EUVI observations of the quiescent prominence around 14:00 UT 
on 2010 December 1 (top row), December 2 (middle row), and December 3 (bottom row). 
(A color version of this figure is available in the online journal.) 
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Fig. 4. — STEREO_B/EUVI observations of the quiescent prominence around 14:00 UT 
on 2010 December 4 (top row), December 5(middle row), and December 6 (bottom row). 
The white and black contours refer to the positive and negative polarities as observed by 
SDO/HMI at 14:00 UT on 2010 December 6. (A color version of this figure is available in 
the online journal.) STEREO_B observations of the prominence at 193 A and 304 A within 
24 hours before the eruption are also available (video 2) in the electronic edition of the 
Astrophysical Journal. 
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Fig. 5. — KSO Ha observations of the quiescent prominence from December 2 to December 
6. 
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Fig. 6.— SDO/AIA (left column, 193 A) and STERECLB/EUVI (right column, 195 A) 
observations of the dense column structure on 2010 December 6 prior to the eruption. (A 
color version of this figure is available in the online journal.) 
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Fig. 7. — SDO/AIA 171 A observations of loops surrounding the prominence prior to the 
eruption. Each image is the average of 30 images taken within 6 minutes. All images are 
contrast enhanced using a radial filter technique developed by S. Cranmer (2010, private 
communication). The field of view of each image is 540" x 540". An animation of this figure 
(video 3) is available in the online journal. Note that radial filter technique is not applied 
for the animation. 
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Fig. 8. — Multi-channel observations of the cavity structure around 18:09 UT on 2010 
December 4. The images in the top two rows are provided by SDO/AIA, while Hinode/XRT 
images are shown in the bottom row. (A color version of this figure is available in the online 
journal.) 
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Fig. 9. — SDO/AIA (left column) of the post-eruption loops at 335 A on 2010 December 7. 
(A color version of this figure is available in the online journal.) 
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Fig. 10.— SDO/AIA (left column) and STERECLB/EUVI (right column) observations of 
the prominence after the eruption. The images at the top, middle, and bottom rows are 
taken at 304 A, 171 A, and 193/195 A on 2010 December 7, respectively. (A color version 
of this figure is available in the online journal.) 
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Fig. 11. — (a) Longitude-latitude map of the radial component of magnetic field in pho- 
tosphere in the HIRES region of the model. The zero-point of the longitude corresponds 
to the central meridian on 2010 December 10 at 14:00 UT. The blue curve shows the path 
along which the flux rope is inserted into the model. Extra axial fluxes are added to (right) 
or removing from (left) the flux rope along places indicated by the yellow bars, (b) The 
magnetic map shown in (a) (red and green contours) is overlaid on the HMI magnetogram 
taken at 14:00 UT on 2012 December 6 (black and white image). The path along which the 
flux rope is inserted is shown as red curve in this image. Note that there is no systematic 
change in the flux distribution from December 6 to December 10. 
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Fig. 12. — Three-dimensional magnetic models of a quiescent prominence observed with 
SDO/AIA on 2010 December 6 at 14:00 UT. The top panels show results from two NLFFF 
models constructed using the flux rope insertion method. For the model in panel (a) the 
axial flux of the flux rope $ ax j = 2 x 10 20 Mx, and the initial poloidal flux F po i = 1 x 10 10 
Mx cm -1 , which produces a relatively small degree of twist. Panel (b) shows a model with 
the same axial flux but larger twist (F po \ = 2 x 10 10 Mx cm -1 ). The colored curves are 
magnetic field lines in and near the flux rope, the blue features indicate field-line dips, and 
the red/green contours show the photospheric flux distribution (B r = ±71 G, respectively). 
The red arc is the east solar limb, and the black lines are the latitudinal boundaries of the 
HIRES domain. These features are projected onto the plane of the sky as seen from the 
SDO spacecraft at 14:00 UT. Panels (c) and (d) show AIA images taken in the 193 A and 
304 A passbands. The path along which the flux rope is inserted is shown as red curve in 
(d). Note that just inside the limb the real prominence is displaced to the south with respect 
to the PIL. This displacement is not reproduced in the NLFFF models. 
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Fig. 13. — Zoomed in view of the field line dips from Model 1 as shown in Figure [T2a near the 
east limb (red arc). The red and green contours refer to the photospheric flux distribution. 
The red curve refers to the path along which the flux rope is inserted. 
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Fig. 14. — Comparison of the model with SDO/AIA and STEREO_BB/EUVI observations of 
the coronal arcade surrounding the prominence cavity (model with poloidal flux F po i = 10 10 
Mx cm -1 ). The top panels show the same set of field lines (colored curves) seen from two 
different viewpoints: (a) SDO and (b) STERECLB on 2010 December 6 at 14:00 UT. The 
red and green contours show the photospheric flux distribution (B r = ±71 G, respectively), 
and the red arcs indicate the east and west limbs. The bottom panels show images in (c) 
AIA 193 A and (d) EUVI 195 A, which are dominated by Fe XII 195 A emission. The 
images show thin threads and fans that presumably are aligned with the local magnetic 
field. On the quiet Sun the emissivity in the Fe XII line drops off rapidly with height (due 
to gravitational stratification of the plasma), so only the plasma on the lower parts of the 
field lines is observable. Therefore, we overplot the lower parts of the field lines that lie at 
heights h < 0.07 R (blue line segments). Note that the orientation of these line segments 
is consistent with the observed fine structures, indicating the model reproduces the fields in 
the coronal arcade. 




Fig. 15. — Magnetic field strength at dips in the field lines from Model 1 are plotted as 
function of the height of the dips above the photosphere. The prominence plasma is assumed 
to be located as such dips. All positions along the filament channel are included in this plot, 
except the two ends of the flux rope where the rope is anchored in the photosphere. Different 
colors represent the field strength at different longitudes as shown in Figure ITTa . 
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Fig. 16. — Top view of the field line dips (blue features) from Model 1 and Model 3. Model 
3 has the same parameters as Model 1 except the filament path (red curve). The red and 
green contours refer to the photospheric flux distribution. 




Fig. 17. — Comparison of the observed cavity and magnetic field lines from Model 1. (a) 
SDO/AIA image of the cavity at 18:08 UT on 2010 December 4. The red and green contours 
refer to photospheric magnetic fields observed by SDO/HMI. (b) The same image as in (a) 
overlaid with magnetic field lines (color lines) surrounding the flux rope from Model 1. 



